Competing supersolids of ultracold Bose-Bose mixtures in a triangular optical lattice 
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We study the ground state properties of a frustrated two-species mixture of bosons in an optical 
triangular lattice, as a function of tunable amplitudes for tunnelling and interactions. By combining 
three different methods, a self-consistent cluster mean-field, exact diagonalizations and effective 
theories, we unravel a very rich and complex phase diagram. More specifically, we discuss the 
existence of three original mixture supersolids: (i) a commensurate with frozen densities and spin- 
like color supersolidity, in a regime of strong interspecies interactions; and (ii) when this interaction 
is weaker, two mutually competing incommensurate supersolids. Finally, we show how these phases 
can be stabilized by the quantum melting of peculiar insulating parent Mott states. 



Ultracold atoms on optical lattices form one of the 
most fascinating fields of research in nowadays condensed 
matter physics. These systems are highly versatile, 
thanks to the possibilities to control the lattice dimen- 
sion and geometry, the statistics of loaded atoms, their 
mutual interactions, as well as to create atom mixtures. 
Therefore, an incredible amount of exotic phases has been 
reported in the literature these last years, both experi- 
mental and theoretical, such as unconventional Mott in- 
sulators, superfluids [l(, Bose metals [2|-|4|| or supersolid 
phases. In this last case, the system enters a phase com- 
bining both crystalline order and superfluidity, and typ- 
ically arising from the quantum melting of a Mott insu- 
lator; this phenomenon is at the origin of intense scien- 
tific activities and debates. Experimentally, supersolidity 
is expected to be realized in solid 4 He [5|, but neces- 
sary conditions for continuous symmetry breaking led to 
question this observation |fj|, which was eventually con- 
tradicted by more recent experiments [7(. However, su- 
persolids are easier achieved on lattice systems; whereas 
on a square geometry with nearest neighbor interactions 
a soft-core description is required [8[, they can also be 
found in hard-core bosonic models when frustration is in- 
duced by cither further neighbor interactions [9[ and/or 
lattice geometry, e.g. triangular |ld-|l2l]. In this context, 
the recent progress in optical lattice engineering gives 
good hopes for observing supersolidity on non-bipartite 
lattices [13|, |l4|. One of the most promising directions 
is the study of mixtures, in particular those with several 
species of bosons, either heteronuclear as 41 K- 87 Rb [L5[ 
or homonuclear with two different spin states, e.g. |1, — 1) 
and |2, — 2) states of 87 Rb [lg|. Another way to realize 
mixtures is to consider, as we will do in this paper, bi- 
layer s yste ms with interactions but no hopping between 
layers [1 71 ] - Such systems allow for an even broader va- 
riety (compared to single-species counterparts) of quan- 
tum phases, depending on the intra- and inter-specics 
interactions, and on the dimension and lattice connec- 
tivity. Compared to an already rich literature on mix- 
tures in ID (e.g. [l8j]) only few works have focused on 
2D cases, for instance on square [19J] and triangular [2C 



22j lattices. In all those systems, contact and dipolar 



interactions 17|, |23|, |2J] can be tuned using Feshbach res- 
onances and by an external electric field [2 51 ] respectively 
Since both interactions, considered independently, tend 
to favor different types of Mott insulators, e.g. density- 
homogeneous or density wave [261 ] . the competition be- 
tween both can, along with quantum fluctuations, trigger 
various unconventional phases 
frustrated lattices. 
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FIG. 1: Left: Triangular optical bilayer lattice as can be real- 
ized by a set of laser beams, and hosting two types of bosons 
labelled a (blue) and b (brown). Right: Sketch of the two- 
component supersolids with 3-sublattice structure: (a) the 
n — 1 spin-like color supersolid (CSS), (b) the 3-supersolid 
(3SS); and (c) the bosonic pinball (BPB), all described in the 
text. The larger the disc, the higher the density. 



In this Letter, we study bosons loaded in an opti- 
cal bilayer triangular lattice, or cquivalently a two-color 
(species) bosonic mixture on a monolayer (see Fig. Q]). 
The richness of the phase diagram of the bosonic ex- 
tended Hubbard model on the triangular lattice, specifi- 
cally with hard core constraint [12(, allows us to expect 
even more exotic physics when the inter-color interac- 
tions are switched on. Indeed, these can mimic efficiently 
the frustrating effects of dipolar couplings that can arise 
in cold atom systems [id . Il2l . 1 17 1 , whereas further range 



interactions have minor effects in the corresponding den- 
sity regimes. Here, we address the following issues: (i) 
how lattice frustration impacts on the stabilization of 
non-conventional phases of such mixture and (ii) how 
on-site and nearest neighbor interactions compete. To 
achieve these goals, we use cluster mean field theory 
(CMFT) and exact diagonalizations (ED) on periodic 
clusters, as well as perturbative approaches. The for- 
mer method has been employed successfully in many one- 
component bosonic systems on various lattices such as 
triangular [12J l26l |31| , Cairo-pentagonal [27| and hexago- 
nal [28j . It is thus expected to be also very efficient in the 
present two-component (bilaycr) model. Our main find- 
ings are the presence of three original two-component su- 
pcrsolid (SS) phases: two incommensurate with density 
modulations, and one commensurate with frozen densi- 
ties and spin-like color supersolidity. 

Model and method- We study a two-component ex- 
tended Bosc-Hubbard model on a triangular lattice: 



lation function S a b = \({n a — ni,)(—k)(n a — nb)(k))\ 2 . We 



used the Fourier transform n a (k) 
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equivalent to a one-component bosons on a bilayer tri- 
angular lattice since no inter-site/inter-color couplings is 
considered. H a is the one-color Hamiltonian (a = a/6), 
with fi a the chemical potential, t a and V a respectively 
the nearest neighbor hopping and interaction amplitudes. 
For consistency, an extra term U aa Ej n ia( n ia ~ l)/2 
describing the on-site repulsion between bosons of same 
color (soft core) has to be included in H a . However, we 
are primarily interested in the effects of frustration in- 
duced by inter color couplings U. These are typically 
reached in the strong interaction limit where the domi- 
nant coupling U aa imposes hard core constaints ni a < 1 
. All the phases reported in this paper can be stabilized 
for U aa > 3V; from now on, b ia is the creation operator 
of a hard core boson of color a at site i. In addition, by 
sampling the parameter space, we found that the most 
original two-component phases arise when the system is 
a/6-symmetric (t a = t, V a = V and // Q = /i) [29j. This 
case, considered in this work, allows by the particle-hole 
canonical transformation b' ia = b ia the use of a mapping 
between cases fj,* > and ^* < 0, with fi* = /u— 3V— U/2 
a rescaled chemical potential. 

We compute the ground states of H a b on clusters up 
to TV = 12 sites by either ED or CMFT methods. In 
the latter case, the Ni (N e ) internal (external) bonds 
are treated exactly (at the mean field level). See Refs. 
12L [30, |3l| for more details. The mean-field parameters 
determined self-consistently are the densities n, a = {rn a ) 
and the superfluid fractions (SF) <j>i a = (bi a ). In addi- 
tion, to evidence 3-fold symmetry breaking, we define the 
order parameter M a = |n Q (A:)| and the color-color corre- 



point k = (47r/3,0), corner of the Brillouin zone. As in 
Ref. [31|, we define the scaling parameter A = N%/(3N) 
which quantifies finite boundary effects. Hence, the Ther- 
modynamic Limit (TL) is achieved for A — > 1 (infinite 

lattice). 
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FIG. 2: Phase diagrams as obtained from CMFT calculations 
on a iV = 3 cluster at t = 0.15 and V — 1. Top: Com- 
plete range of fj,* showing two-component superfluids, 2/3 and 
4/3 plateaus, XXZ physics domain and pure two-component 
physics not reminiscent to the phases of the one-specie boson 
system (white rectangle) . Bottom: zoom of the white rectan- 
gle and all new two-component phases (see text for acronyms) . 
Solid (dashed) lines are phase boundaries subsisting (vanish- 
ing) in the Thermodynamic Limit (TL). 



Overview of the phase diagram- Let us first focus 
on the (fi* , U) phase diagram in which novel two-color 
phases - detailed along the paper - emerge, illustrated 
here for t/V = 0.15. For these parameters, the one-color 
Hamiltonian H a is known to present a rich phase diagram 
with either empty /full, homogeneous superfluid, V3x v3 
solid at density 1/3 or supersolid phases [lfj. When U 
is switched on, the correlations between the two colors 
increase, and resulting GS can be either a product of 
two one-color phases or color-entangled states as depicted 
in the complete phase diagram shown in Fig. [5] (i) As 
/j* < increases, the first non-trivial phase encountered 
is a homogeneous two-color superfluid with <p a = (\>\, ^ 
(4>a = j? J2i <£«*)■ When n a and n b are high enough, V 
and U terms become costful and the system enters a 3- 
fold ordered insulator characterized by <f> a = (f>b = and 
a total density n = 2/3 (n a = ni, = 1/3), the two-color 
counterpart of the \/3 x \/3 phase, (ii) At even larger 
density, various supersolids with 3-sublattice structures 
are stabilized. Finally, a phase with n = 1 is stabilized 
for U > 1.2V. For fi* > 0, we obtain an equivalent phase 
diagram thanks to the particle-hole symmetry. The most 
original of these two-color phases are detailed hereafter. 
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FIG. 3: Strong U regime characterized using either N — 12 
(ED: a,d) and N = 6 (CMFT: b,c) clusters, (a) n and SU as 
a function of ^*/V for U = 2V and t = 0.15V. The n = 1 
plateau and the finite S a b evidence the CSS; (b) n and (/> = 
0o + <t>b as a function of [i* /U distinguishing the CSF (V = 0) 
and the CSS (V — 0.2(7) regimes; (c) quantum melting of 
the CSS to the BPB phase (see text); (d) comparison of the 
kinetic energy gain e* ~ —kt 2 /U for H a b (filled) and -Hxxz 
(empty) for three (£, V) sets. 



only off-diagonal order [32], [33j. In the present context, 
these phases correspond respectively to a color supersolid 
(CSS) and a color superfluid (CSF). While the density 
is uniform, the CSS has a (2m z , — m z , — m z ) structure 
where m z quantifies the color disproportion on each sub- 
lattice. As depicted in Fig. [3J for \/j,*\ < V, both ED 
and CMFT are consistent with the CSS scenario; indeed 
the finite structure factor S a b and M = (M a + Mb)/ 2 for 
t/U < 0.17 signal long-range correlations and the corre- 
sponding XXZ couplings verify J z /J± = 1 + UV/2t 2 > 
20.5. In contrast, for V = and t/U < {t/U) c ~ 0.05(1) 
we find in vicinity of \x* =0 the spatially uniform CSF 
predicted in the XXZ model, as illustrated in Fig. EJb). 
Finally, Fig. (3Jd) shows that the kinetic energy gain e* ~ 
—kt 2 /U w.r.t. the electrostatic contribution N(V — (J,) is 
well reproduced by the XXZ model in both regimes; this 
validates our approach. 
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FIG. 4: Color-resolved observables {n a ,M a ,<f> a ). Each color 
corresponds to either continuous or dashed line. All data 
come from a N = 6 CMFT with fixed t/V = 0.15 and U = 
0.5V (left) or 2V (right). The arrow points out a tiny region 
which disappears in the TL (see text). The shaded regions 
correspond to the different phases of Fig. [l] 



Commensurate spin-like phases- For strong inter-color 
repulsion U 3> V, t, |/x*| (triangular uppermost domain 
on Fig. [H-up) the system is Mott insulating, with <f> a = 
and n = 1. Indeed, U imposes the local constraint of 
single occupancy defined as n a + nb = 1. This is reflected 
in ED by the n = 1 plateau on Fig. ^,b) and in CMFT 
by the n a = 1/2 plateau on Fig. Upright). In order to 
better describe this regime, we define spin 1/2 operators 
erf = (rii a — riib)/2 and oy = b ia bn,. At second order 
of the perturbation theory, we obtain an effective XXZ 
model 



H 



XXZ 
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where J± = At 2 /U and J z = 2V + At 2 /U . Interest- 
ingly, this model predicts a (2m z , — m Zl — m z ) SS with 
diagonal (er 2 ) ^ and off-diagonal (er + ) ^ order pa- 
rameters when J z /J± > 4.6(1); otherwise, a SF with 



Incommensurate supersolids- In contrast to the above 
mentioned CSS in which supersolidity comes from the 
color spin-like degrees of freedom, in SS phases discussed 
here n a may be incommensurate and they have finite </> a . 
From the color-resolved observables shown in Fig. 01 we 
identify three such SS phases, (i) For large U/V = 2 
(right column) and fi*/V about ±1.0(2), the finite M a 
indicates a 3-fold order (3FO) for both colors, while only 
one component has a non-zero <p a . This phase, dubbed 
SS/3FO in Fig. [21 is a rearrangement of two one-color 
phases of H a , the -\/3 x -\/3 solid and the SS obtained 
from its melting [lfj, [l2| . Two sublattices are (almost) 
filled by a and b bosons respectively, while on the remain- 
ing, an incommensurate density for bosons of one color 
accounts for superfluidity. Within this structure, the re- 
pulsion energy oc U is minimized thanks to the localiza- 
tion of bosons of a single color, (ii) For small U/V = 0.5 
(left column) and [i* /V about ±0.5(3), the phase has a 
finite 4> a for both colors as well as a 3-fold symmetry 



breaking. This phase, called 3-fold SS (3SS) is the first 
example of a collective two-color supersolid. As shown 
in Fig. HJb), it can be seen as a superposition of two 
one-color SS (these would be obtained at U = 0); both 
colors contribute symmetrically to superfluidity. Note 
that this 3SS is obtained from the parent n = 2/3 solid 
by a defect condensation upon doping as /i* increases. As 
in this parent state, the weak inter-color coupling merely 
forces the localized a and b bosons to occupy distinct sub- 
lattices, (iii) The most remarkable incommensurate SS 
phase is achieved by the quantum melting of the density- 
uniform n = 1 CSS. This original two-color SS is called 
the bosonic pinball (BPB) due to a structure very similar 
to its fermionic counterpart with similar interactions, the 
pinball liquid. It has almost one localized electron per site 
on one sublattice (pins) and a metallic behavior on the 
remaining hexagonal lattice (balls) |26j, |34|, |35[ . Here, the 
particles are bosonic and the colors play the role of the 
spins. The structure is depicted on FigQJc). One sub- 
lattice, forming a triangular super-lattice with fii a + nu, 
are close to 1, is filled in majority by one type of bosons. 
The two-color supcrfluid character is carried by the re- 
maining bosons on the complementary hexagonal lattice. 
This is shown in the BPB region on Fig. |4|left column) 
and Fig. El by a coexistence of a solid M a -£ and su- 
pcrfluid (p a 7^ orders in both colors. This lattice sym- 
metry breaking is reminiscent of the parent CSS phase, 
the quantum melting of which involves a condensation of 
two types of defects coming from doubly-occupied and 
empty sites. Finally, unlike the 3SS and SS/3FO, the 
BPB can be stabilized at the n = 1 commensurability 
when (i* = 0. In Fig. [5] are displayed the cluster depen- 
dency of M and <f> for two examples illustrating the BPB 
and the 3SS which confirms their existence at the TL as 
A — >• 1 (inset). 




FIG. 5: M and (j) as a, function of fj,* for t - 0.15V and U = 
0.5V obtained by CMFT on clusters of N — 3 (continuous), 
6 (dashed) and 12 (symbols) sites. (Inset) cluster scaling as 
function of A (see text) of M (filled symbols) and <jf>/2 (empty 
symbols), for the BPB at \fi*\ = t (squares) and the 3SS at 
|/i*| — 3t (circles). 



Phase transitions- We checked that the phase bound- 
aries are only weakly affected by size effects. The phase 



diagram of Fig. [2] is then representative of the TL. It is 
possible to extract an estimate of certain of these tran- 
sitions when defects proliferate as their energy cost van- 
ish (condensation picture), in the strong correlation limit 
U, V > t [3|, as it is the case for the 2/3 -> 3SS tran- 
sition. Consider a defect consisting of an extra boson 
(say of color a) inserted in the empty sublattice of the 
2/3 crystal with the energy cost 3V. Via second or- 
der processes, such a defect can hop to second neigh- 
bors with an amplitude -t e g = -t 2 /U - t 2 /V J33j- In 
this limit, the defects condense and lead to supersolidity 
for n c = W - 6t off , e.g. ^ C /V ~ -0.65 for U = 0.5V 
in good agreement with the results of Fig. [U This de- 
fect condensation mechanism implies a gauge symmetry 
breaking additionally to the already broken lattice sym- 
metry of the 2/3 phase. This indicates that the tran- 
sition is of second order, as confirmed by the absence 
of discontinuities in <j> and M as function of /i (Fig. [5]). 
This corresponds to the standard picture of defect con- 
densation accounting for an incommensurate SS; the SF 
density (here inhomogeneous) is carried by defects on top 
of density modulations [lCJ. In contrast, the transitions 
between (i) homogeneous SF and 2/3 solid, and (ii) the 
BPB and 3SS are found of first order, characterized by 
hysteresis in the CMFT [38( . In case (ii) for example, 3SS 
and BPB have distinct broken symmetries as displayed 
inFig.m:b,c). 

Conclusions- We study interacting hard-core bosons 
in an optical bilayer triangular lattice equivalent to a 
two-color bosonic mixture on a monolayer, by combining 
both ED and CMFT methods. We focus on a region of 
parameter space in which novel phases arise due to the 
competition between frustration and quantum fluctua- 
tions. Within a very rich and complex phase diagram, 
we evidence three original mixture supcrsolids with 3- 
site unit cell structures: (i) a spin-like commensurate 
phase termed the color supersolid (CSS), characterized 
by an effective XXZ model with frozen densities and su- 
persolidity carried by the color degrees of freedom, (ii) 
two mutually competing incommensurate phases arising 
from the quantum melting of parent states, either the 
2/3-solid or the original CSS, via defect condensation 
mechanisms. Melting the 2/3 phase leads to a 3-fold 
supersolid (3SS) with three inequivalcnt sites. The most 
interesting supersolid comes from the melting of the CSS 
and directly results from strong effects of inter-color cou- 
pling. The inner structure of this bosonic pinball (BPB) 
is very reminiscent to that of its fermionic counterpart. 
These striking features found in the hard-core limit and 
by truncating long-range to n.n. interactions tunable ex- 
perimentally, are expectable in a broader variety of situa- 
tions (dipolar interactions, sot-core bosons) which should 
stimulate further investigations. 
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